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We investigate potentials between pairs of static-light mesons in Nf = 2 Lattice QCD, in different 
spin channels. The question of attraction and repulsion is particularly interesting with respect to 
the X(3872) charmonium state and charged candidates such as the Z + (4430). We employ the 
nonperturbatively improved Sheikholeslami-Wohlert fermion and the Wilson gauge actions at a 
lattice spacing a « 0.084 fm and a pseudoscalar mass mps ~ 760 MeV. We use stochastic all-to-all 
propagator techniques, improved by a hopping parameter expansion. The analysis is based on the 
variational method, utilizing various source and sink interpolators. 
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1. Introduction 

We numerically determine ground and excited states of static-light mesons (3$ = Qq) as well as 
intermeson potentials between pairs of static-light mesons, 3§{r)S8{Qi) and ^(r)^(0), with a static 
quark-quark (or quark-antiquark) separation r = |r| = Ra. a denotes the lattice spacing, Q a static 
colour source and the positions of the light quarks q are not fixed. For large heavy quark masses the 
spectra of heavy-light mesons are determined by excitations of the light quark and gluonic degrees 
of freedom. In particular, the vector-pseudoscalar splitting vanishes and the static-light meson 38 
can be interpreted as either a B, a B* , a D or a D* heavy-light meson. 

Static-light intermeson potentials were first evaluated on the lattice by Michael and Pennanen 
in the quenched approximation [1] and with Nf = 2 Sheikholeslami-Wohlert sea quarks [2]. A more 
detailed quenched study can be found in ref. [3] and the un-quenched case is revisited with twisted 
mass fermions by Wagner in these proceedings [4]. Meson-antimeson potentials were computed 
in ref. [2] and, with Wilson sea quarks, in ref. [5]. Graphically, the quark line diagrams that we 
evaluate can be depicted as, 



0- Of M 



(i.i) 
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where the static-light correlation function is shown on the left, and those for the meson-meson 
and meson-antimeson cases in the centre and on the right, respectively. Straight lines represent 
static quark propagators and wiggly lines light quark propagators. In this paper we neglect explicit 
meson exchange (i.e. "box" and "cross") diagrams. The analyses of these as well as of a larger 
lattice volume are in progress. This means that here we only consider the isospin 1 = and 
the 1=1 £$S$ combinations. 

The static-light correlation function in Euclidean time t = Ta is given by, 



C(t) = (0\(Q a G a p q p ) x+d (q r G yS Q 5 ),|0) = ( Tr 



1 1 x+kcAA O x+t4,x 
k=0 




The trace is over colour (not displayed) and Dirac indices, (-) v indicates the expectation value 
over gauge configurations and /t denotes a unit vector in fi -direction. We average over all possible 
source points x/a € {1, .. . ,L a } 3 x {1, . . . ,L T } to reduce statistical errors where x = (x,^). The 
correlator is automatically zero-momentum projected since x is the same at source and sink, due 
to the static propagator. MZ. = {q y q x ) is the propagator for the light quark q on a given gauge 
configuration and U xll is the gauge link connecting the lattice site x with x + ajx. The absence 
of the spin in the static propagator necessitates the (1 + /4)/2 Dirac projection of the (fermionic) 
static-light "meson" to fix the parity P. This is very similar to baryonic correlation functions where 
a spin j source is created by three (rather than one) light quarks. Meson-(anti)meson correlation 
functions can be obtained by combining the above correlator with another one that is spatially 
shifted by a distance r, before taking the gauge average. 
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2. Representations and classification of states 

In the continuum limit, the static-light states can be classified according to fermionic repre- 
sentations J p of the rotation group 0(3). At vanishing distance r = the and states can 
be characterized by integer J p and J PC quantum numbers, respectively. However at r = |r| > 
the 0(3) (or 0(3) (g)^ 7 ) symmetry is broken down to its cylindrical Dooh subgroup. The irreducible 
representations of this are conventionally labelled by the spin along the axis A, where £,11, A refer 
to A = 0, 1 , 2, respectively, with a subscript f] = g for gerade (even) PC = + or T] = u for ungerade 
(odd) PC = — transformation properties with respect to the midpoint. All A > 1 representations are 
two-dimensional. The one-dimensional £ representations carry an additional o v = ± superscript 
for their reflection symmetry with respect to a plane that includes the two endpoints. 

To create states of different J p ( c ) we use operators & that contain combinations of Dirac y- 
matrices and covariant lattice derivatives V[U] that act on a fermion spinor q as, 

V^.v = U Xjfl q y+a p, - U x ^q x - a (i , where U x = U*_ afL ^ . (2. 1) 

On the lattice the continuum rotational symmetry is broken and the groups 0(3) and Dooh need to 
be replaced by their finite dimensional subgroups Oh and D 4h , respectively. We label fermionic 
representations of the octahedral group Oh as Oh'- For fermionic representations of Do»h that we do 
not need in the present context, see ref. [6]. It is well known, see e.g. ref. [7], that the assignment 
of a continuum spin to a lattice result can be ambiguous, in particular for radial excitations because 
a given Oh representation can be subduced from several continuum 7s. For instance, 

GK- J H^J=^, 5 -,..., Ai<-/ = 0,4,-, 7\ ^7= 1,3,4,.... (2.2) 

For A < 2 the mapping of continuum Dooh onto discrete D4h representations is more straight for- 
ward. Hence in this case we adopt the continuum notation only. 

The operators that we used to create the static-light mesons are displayed in table 1 (see, e.g., 
ref. [8]). The intermeson potentials were obtained by combining two static-light mesons of differ- 
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Table 1: Operators and representations for static-light mesons. In the last column we display the / for a 
heavy-light meson, obtained by substituting the (spinless) static source by a heavy fermion. 

ent (or the same) quantum numbers. This can be projected into an irreducible Dooh representation, 
either by coupling the light quarks together in spinor space [4] or by projecting the static-light 
meson spins into the direction f of the static source distance, by applying g(l ± iy^y ■ f), and 
taking appropriate symmetric (A s = 1) or antisymmetric (A s = 0) spin combinations. These two 
approaches can be related to each other via a Fierz transformation. For the preliminary results pre- 
sented here we have not yet performed this projection and different representations will mix. The 
analyzed operators and the corresponding representations are listed in table 2. 
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Table 2: Operators and continuum representations for the meson-meson (SSSS) and meson-antimeson (9838) 
potentials. 



For / > and r > the irreducible representations of 0(3) split up into two or more irreducible 
representations of Dcoh- For instance the angular momentum of the P + operator within our 1~,2~ 
r = state can be perpendicular or parallel to the intermeson axis. For the axis pointing into the 3- 
direction, we call the S&P+ operator 75 <g) (yiVi — 72 V2) "parallel" (||) and the other combinations 
"perpendicular" (_L). The _L state has no angular momentum pointing into the direction of the axis 
and hence only couples with the light quarks to E+ and IT,,. Vice versa, the || operator can only 
create H u and A u states but not the £+. 

3. Simulation and analysis 



volume x L t 
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Table 3: Lattice parameters. 



We employ Nf = 2 Sheikholeslami-Wohlert configurations generated by the QCDSF Collab- 
oration [9]. The parameter values are listed in table 3, where the scale is set using tq(P,k) = 
0.467 fm. The pseudoscalar mass corresponds to its infinite volume value. We use the Chroma 
software system [10]. 

To achieve high statistics in the evaluation of the diagrams eq. (1.1), all- to-all propagators 
need to be computed. This is done using stochastic estimator techniques, see ref. [11] and ref- 
erences therein. We generate 300 complex Z2 noise sources and apply the hopping parameter 
expansion to reduce the stochastic variance [12, 5, 11]. Furthermore we enhance the signal over 
noise ratio by employing a static action with reduced self-energy [5]. This is done by applying 
one stout smearing step [13] with the parameter p = 1/6 to the temporal links, used to calculate 
the static propagators. Wuppertal smearing [14] is applied to the source and sink operators, where 
we employ spatially smeared parallel transporters [5] with the parameters n lier = 15, a = 2.5. The 
Wuppertal smearing hopping parameter value fc w = 0.3 is combined with three iteration numbers 
Mter £ {16,50, 100}. Masses are then extracted from the resulting correlation matrices, by means 
of the variational method [15], solving a generalized eigenvalue problem. Errors are calculated 
using the jackknife method. 
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4. Results 

The eigenvalues A ^(f,fo) of the generalized eigenvalue problem [15], 

Cij{t)uf =^(t,to)Qj(to)uf , 

are fitted to one- and two-exponential ansatze, to obtain the kth mass. The appropriate values of to 
and the fit ranges in t are determined from monitoring the effective masses, 



% a {t + a/2)=a-'\n 



A(*)(f + a,/d) 



On the left hand side of figure 1 , we display the effective ground state energy levels E e s for to = 2a 
of the 75 (g> 75 SSSS system for different distances R = r/a. In the limit r — > oo these will approach 
the sum of two 4 static-light meson effective masses (dotted curve). At short distances we see 
attraction in this channel. The quality of the effective mass plateaus deteriorates with decreasing 
distance since the wavefunction of such interacting states becomes more than the mere product 
of our two static-light meson interpolators. We should also keep in mind that so far we did not 
perform the singlet spin projection and hence there will be additional pollution from the Yl g state, 
see table 2. For the example displayed, we perform two-exponential fits to the t /a £ {4, ... , 10} 
data to obtain the masses. 
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Figure 1: Effective groundstate masses of the 75 £g> 75 operator at different distances (left hand side). The 
dotted line corresponds to twice this mass for a single ^ + static-light meson. The intermeson potential V(r) 
for the combinations 75 <£> 75 (E+ £+') and 75 ® 1 (E^") (right hand side). 



We define intermeson potentials as the difference between the meson-meson energy levels and 
the r — > 00 two static-light meson limiting cases: 



%,^ 2 ( r ) = E &i3h ( r ) ~ ( m m + m a 2 ) 



0. 



(4.1) 



The results for the groundstate (£j) and the first excited state (£j ) of the 75 ® 75 operator and the 
Z u groundstate of 75 (8) 1 are plotted on the right hand side of figure 1. In all these 1 = channels 
there is attraction of the order of 50 MeV at a distance of 0.2 fm. 

In figure 2 we display the ground state and first excited state energy levels for the SSSS meson- 
antimeson case in the 75 (g) 75 channel. The two horizontal lines correspond to twice the ground state 
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Figure 2: Effective masses for the 75 <g> 75 38 SB ground and first excited states (red symbols). The horizontal 
lines denote sums of two static-light meson masses while the black curve and symbols are the static QQ 
potential, offset by the scalar «o meson mass. 



mass of the ^ + static-light meson and to the sum of its ground and its first excited state masses, 
the expected r — > 00 limits. At first sight there appear to be very substantial short distance attractive 
forces in this channel. However, states consisting of a QQ static potential and a scalar 1 = 1 particle 
will have the same quantum numbers. For our lattice parameters the ao meson is the lowest such 
state, with masses of two pseudoscalars as well as of a P-wave vector lying higher. We include the 
sum of these two masses in the figure. The ground state 333§ energy still lies below this level but 
to decide whether we effectively see the sum of ao and the static potential and to disentangle which 
/ = 1 energy level is the lowest one, interactions of the ao with the static potential will have to 
be taken into account. We hope that simulations on a larger volume will help to clarify this. 



5. Conclusions 



We investigated interactions between pairs of static-light mesons and found attraction in the 
I = S = sector. Meson-antimeson potentials are also very interesting with respect to charmonium 
threshold states [16] (DD molecules or tetraquarks) but difficult to disentangle from mesons that are 
bound to static-static states (hadro-quarkonium [17]). We are in the process to extend our study to 
1=1 and to / = 3$ S3 states and to disentangle = ±1 from 5f = states. Also simulations 
on a larger volume are in progress. 
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